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Abstract 

We establish that the static height fluctuations of a particular growth model, the PNG droplet, 
converges upon proper rescaling to a hmit process, which we call the Airy process A{y). The 
Airy process is stationary, it has continuous sample paths, its single "time" (flxed y) distribution 
is the Tracy- Widom distribution of the largest eigenvalue of a GUE random matrix, and the 
Airy process has a slow decay of correlations as y~'^. Roughly the Airy process describes the 
last line of Dyson's Brownian motion model for random matrices. Our construction uses a 
multi-layer version of the PNG model, which can be analyzed through fermionic techniques. 
Specializing our result to a flxed value of y, one reobtains the celebrated result of Baik, Deift, 
and Johansson on the length of the longest increasing subsequence of a random permutation. 



1 The PNG droplet 



The polynuclear growth (PNG) model is a simphfied model for layer by layer growth [|I|, §]. 
Initially one has a perfectly flat crystal in contact with its supersaturated vapor. Once in 
a while a supercritical seed is formed, which then spreads laterally by further attachment of 
particles at its perimeter sites. Such islands coalesce if they are in the same layer and further 
islands may be nucleated upon already existing ones. The PNG model ignores the lateral lattice 
structure and assumes that the islands are circular and spread at constant speed. 

In this paper we study the one-dimensional version in the particular geometry where nucle- 
ation only above the ground layer [—t, t] is allowed. To be precise: at time t the height is given 
by the (random) function h{x,t), x G M. One requires h{x,t) = for |a;| > t, in particular 
h{x, 0) = for X G M. At a given time t, h{x, t) is piecewise constant, takes nonnegative integer 
values, and has jumps of size ±1 only. The jumps are called steps and we distinguish between 
up-steps (jump size 1) and down-steps (jump size —1). The dynamics has a deterministic piece, 
according to which down-steps move with velocity 1 and up-steps with velocity —1. Surface 
steps disappear upon collision. In addition there are nucleation events by which new steps 
are created. Randomly in time, h{x,t) is changed to the new profile h{x,t) such that for the 
increment 6h{x, t) = h{x, t) — h{x, t) one has 5h{x, t) = 1 at some random point x', \x'\ < t, and 
6h{x,t) = otherwise. Immediately after this nucleation event the deterministic evolution is 
followed until the next nucleation. Growing for a while the typical height profile has the shape 
of a droplet, h{x,t) ^ 2\fW^^^ for \x\ < t. Our interest are the statistical properties of the 
deviations from this average shape. 

Warning: For the PNG model one has to specify the step speed and the intensity of the nu- 
cleation events. They can be adjusted to an arbitrary value through a linear scale change of 
space-time (x, t). Geometrically velocity one is distinguished and therefore adopted here. In- 
tensity one for nucleation events seems also natural, but in fact introduces a string of factors 
V^. Therefore we deviate from previous conventions and set the intensity to be equal to 2. 

The one-dimensional PNG model is just one model within the KPZ universality class for 
growth. However, for this model we have very refined statistical information, the most surprising 
breakthrough being the result of Baik, Deift, and Johansson [Q, which states that 

limt-i/3(/^(0,t)-2t) = X2 (1.1) 

f— >oo 

in distribution, cf. for the connection to the PNG model at the end of the introduction. X2 has 
the same distribution as the largest eigenvalue of a x random hermitian matrix (GUE) 
in the limit of A^ ^ oo. As discovered by Tracy and Widom the distribution function 
-^2(2;) = IP(X2 ^ x) is governed by the Painleve II equation. One has F2{x) = e~^^^\ where 
(7" = n^, (^(x) ^ as X ^ 00, and u is the global positive solution of u" = 2u^ + xu (Painleve 
II). Its asymptotics are u{x) ~ —x/2 for x —00, u{x) ~ Ai(x) for x ^ 00 with Ai the Airy 
function. In fact through a simple linear transformation one obtains the height fluctuations 
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for any \y\ < 1, as 

At the next level of precision, one might wonder about joint distributions of h{x,t) at the 
same time, say of {h{xi,t), h{x2,t)}. If Xi = yit and X2 = y2t, yi 7^ y2, then the heights are 
so far apart that they become statistically independent. Thus we have to consider closer by 
reference points. From the KPZ theory one knows that the lateral fluctuations live on the scale 
t^/^. Therefore the natural object is 

y^ry^{h{yt'/\t)-2t)=h,{y) (1.3) 

considered as a stochastic process in y. For large t we have {ht{y)) — — and if in ( p..2| ) we 
replace yt by yt"^^^ we obtain 

lim ht{y) = -y^ + X2 (1.4) 

in distribution, which suggests that ht{y) + y"^ tends to a stationary stochastic process. As our 
main result we establish that this is indeed the case and rather explicitly identify the limit 
process. For reasons which will become clear in the sequel we call the novel limit process the 
Airy process and denote it by Aijj). We refer to Section ^ for its definition and to Section ^ 
for some of its properties. Somewhat compressed, one considers independent fermions in one 
space dimension as governed by the one-particle Hamiltonian 

+ (1.5) 

The fermions are in their ground state at zero chemical potential, which is the quasifree state 
determined by the spectral projection onto {H < 0}. Because of the linearly increasing po- 
tential there is a last fermion, which has the Tracy- Widom X2 as positional distribution. Ex- 
tending to the Euclidean space-time through the propagator e~^^ , the fermions move along 
non- intersecting world lines with some suitable statistical weight. The Airy process A{y) is the 
position of the last fermion at fermionic time y. 

Theorem 1.1 Let A{y) be the stationary Airy process. Then in the sense of weak convergence 
of finite dimensional distributions 

lim ht{y) = A{y) - y\ (1.6) 

t— >oo 



We recover the result in 0] as the special case of the convergence of the distributions for some 
fixed value of y. Even for fixed y we provide an alternative proof based on one-dimensional 
fermionic field theories. 
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As shown in [|], ^] the PNG droplet is isomorphic to a directed polymer with Poissonian 
point potential which in turn is related to the length of the longest increasing subsequence in 
a random permutation. For the sake of completeness we repeat our result in this language. 
We consider Poisson points of intensity 2 in the positive quadrant. A directed polymer is a 
piecewise linear path, uj, from (0, 0) to (t + x, t — x), |x| < t, with each segment of u bordered 
by Poisson points, under the constraint that their slope is in [0, oo]. The length, L{x,t,u!), of 
the directed polymer u is the number of Poisson points visited by u. We set 

L{x,t) = ma.xL{x,t,uj) , (1.7) 

UJ 

where the maximum is taken over all allowed directed paths at a fixed configuration of Poisson 
points and at specified endpoints. Then, in distribution, 

L{x,t) = h{x,t) . (1.8) 

Therefore Theorem p..l| yields the statistical properties of the length of the optimal path in 
dependence on its endpoint at transverse distance yt'^^^ away from {t,t). 

To give a brief outline: In the following section we introduce the multi-layer PNG model, 
whose last layer is the PNG droplet. The multi- layer PNG has the remarkable property that 
the distribution at time t is the uniform distribution on all admissible height lines. Such kind 
of ensemble can be analyzed through Euclidean Fermi fields. Our case maps onto independent 
fermions on the lattice Z with the usual nearest neighbor hopping energy and subject to a 
linearly increasing external potential. The PNG droplet corresponds to the last fermionic world 
line. The convergence of the moments of the multi-layer PNG model in essence reduces to the 
convergence of the discrete Fermi propagator to the continuum Fermi propagator corresponding 
to the Hamiltonian ( |1.5| ). In the final section we establish some properties of the Airy process 
and discuss the two-point function {(h{yt'^/^,t) — /i(0,t))^). 

2 The multi-layer PNG model 

Inspired by the beautiful work of Johansson on the Aztec diamond ^ we enlarge the PNG 
droplet to the multi-layer PNG model. While this looks like a further complication, in fact the 
construction will provide us with a powerful machinery to analyze the statistics of the PNG 
droplet. 

Instead of a single PNG line we now consider a collection of such lines, denoted by h£{x,t), 
£ G N_ = {0, —1, —2, . . . }, taking values in Z, which are arranged in ascending order, h£{x, t) < 
h£^i{x,t), for all X G M, t > 0, cf. Figure |l| with a typical snapshot of the multi-layer PNG 
model at time t. The lines have jump size ±1 and the total number of up- and down-steps 
is assumed to be finite. For each t > one has h£{x,t) = £ for |x| > t. To be definite we 
require he{x, t) — lime\o hi^x ± e, t) G {0, 1}, which means that the height lines are upper semi- 
continuous. The set of all such height line configurations is denoted by A^. The height lines 
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Figure 1: A snapshot of a multi-layer PNG configuration at time t. As a guide to the eye the 
asymptotic droplet shape is indicated. 

evolve in time under a rule which is based on Viennot's geometric construction to prove the 
Robinson-Schensted(-Knuth) correspondence ||^, hence called RSK dynamics by us. The initial 
line configuration is /i^(a:, 0) = £, £ G N_, for all x. Under the RSK dynamics only the region 
\x\ < t is modified. The top line evolves stochastically like the PNG droplet. In the lower 
lying lines the steps move and coalesce according to the PNG rules. However nucleations are 
determined by the annihilation events in the neighboring line above. Thus at time t if in the 
i-th height line a collision of an up-step and a down-step occurs at position x, they disappear 
in this line only to reappear as nucleation at (x, t) in line i — 1. Clearly the RSK dynamics 
respects the ordering hi < /i^+i. 

To describe the configuration space of the multi- layer PNG model at time t > 0, we introduce 
the step positions as coordinates. Let n£ be the number of up-steps in height line i. Since 
hi{—t,t) = h£{t,t), it must equal the number of down-steps in height line i. If ne = 0, 
hii^x, t) = £ for all x G M. If > 0, the position of the j-th up-step in height line £ is denoted 
by yj'~^, —t < yf~^ < ■ ■ ■ < y^^ < t and the position of the j-th down-step in the same height line 
is denoted by yf', —t < yi~ < ■ ■ ■ < y^i" < t. We set n = {uq, n^i, . . . ), and |n| = J2e<o 
For the RSK dynamics the total number of steps, 2|n|, is finite with probability one, since |n| 
equals the total number of nucleation events up to time t, i.e. the number of Poisson points 
in the triangle \x'\ < t', t' < t}. We denote by Tt{n) the set of all step configurations 

{{yf^ 'yf~)i<j<ne) i<,o resulting from an admissible line configuration [hi{x,t))^^^ G A*. Tt{0) 
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is a single point and Tt{n) is naturally embedded in [—t, t]^l"L By definition of A^, if an up-step 
and a down-step in the same height line are at the same location, they represent a nucleation 
and not a collision. Finally Tt = U|n|<oo -'^t(^)- have thus defined the map S : At ^ Tt, 
which we call step map. Clearly, 5" is invertible. The RSK dynamics on Af induces a dynamics 
of steps on Tt- By construction the RSK dynamics stays inside Tt with probability one. 

Remarkably, the distribution of the multi-layer PNG model at time t has a simple structure 
in being the uniform Lebesgue measure on all admissible step configurations. 

Theorem 2.1 Let Wt be the uniform measure on Tt, which means that Wt[Tt{0)^ = 1 and 
Wt \ Tt{n) is the 2\n\- dimensional Lebesgue measure on r((n). Then J wt = Z{t) = exp(2t^) 
and fit = Z{t)~^wt is a probability measure on Tt- If the height lines evolve under the RSK 
dynamics, then fit is the joint distribution of {h£{x,t), x G M, £ G N_} under the step map S. 

Proof. The nucleation events determining the line configuration at time t are a Poisson process 
of intensity 2 in the triangle < t',0 < t' < t}. If {x',t') is a generic Poisson 

point, we label it through the new coordinates {y^,y^) = {x' — {t — t'),x' + (t — t')). A 
Poisson point configuration consisting of points is then given by {y^ ,y~)i=i,...,N, such that 
—t<yf<---<y'^<t and y^ < y^ , i = 1, . . . , N. The set of all such point configurations 
is denoted by At{N), considered subset of We also set = [Jj^^^^At{N), with 

At(0) a single point. inherits from the Poisson process the probability measure i>t, where 
i^t{At{0)) = and z/* \ At{N) = e-^'' dyt ■ --dy^dy^ ■ --dy],. 

Next we define the growth map G : At ^ At. For given {yf iyi)i<i<N £ At{N) we run 
the RSK dynamics to determine the line configuration at time t. Conversely for a given line 
configuration in At we run the RSK dynamics backwards in time which then determines the 
Poisson points corresponding to the nucleation events. Thus G is well defined. If all the line 
configurations with coinciding up- resp. down-step positions in different lines — a set of Lebesgue 
measure zero under the step map — are removed from A^, then G is even bijective. 

From the construction it is obvious that foiSoG: {yf ,yi)i=i,...^N ^ {iyf~^ ■>yf~)i<j<nt) i<,q 
with A^ = |n| one has the following set equalities, {y^, i = 1, . . . , N} = {yf^, j = 1, . . . , n^, £ G 
N_} and {y^ , i = 1, . . . , A^} = {yf', j = ^, ■ ■ ■ ,n£, £ E N_}. Thus the map S o G induces a 
mere relabeling of points. In particular Ut is transformed to fit under S o G. O 

Equipped with Theorem |2]l| the reader may jump ahead to Section ^ where the statisti- 
cal properties of the line ensemble /if are studied. We take a little detour to report on two 
observations of interest. First there is a variant of the multi-layer PNG dynamics which was 
introduced by Gates and Westcott in modeling crystal growth, hence called GW dynamics 
by us. Bulk properties of the GW dynamics are studied in ^ in the context of the anisotropic 
KPZ equation. Gates and Westcott regard the lines hi as contour lines of a crystal surface. 
The crystal is made up of atomic two-dimensional layers stacked along the z-axis. Layer and 
below are completely filled. For layer 1 only the domain {(x, y) G ?/ < ho(x, t)} is filled with 
atoms, and in general, layer —£+1 is filled in the domain {(x,y) G M^l y < hi(x,t)—£}, £ G N_. 
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The crystal is in contact with its supersaturated vapor. If overhangs are not permitted, it is a 
natural modeling assumption that each contour line grows under the PNG rules subject to a 
constraint of no touching. In this sense the GW dynamics is "more stochastic" than the RSK 
dynamics. Despite different rules, the RSK and GW dynamics yield an identical distribution 



for the line configurations at time t. Thus, alternatively, the proof of Theorem could be 
based on the GW dynamics, cf. below. 

A further observation is that the GW dynamics admits a discrete space-time version which, 
as to be discussed at the end of this section, inherits the simplicity of the distribution at time 
t. In fact, the discrete version of the GW dynamics is isomorphic to the shuffling algorithm for 
the Aztec diamond. This provides us with yet another approach to Theorem |2.1| , namely to 
take the continuum limit of its discrete analogue. 

Let us start with the GW dynamics. The top line evolves as the PNG droplet. With 
nucleation rate 2 pairs of up- and down-steps are generated in the forward light cone of the origin 
and move apart with velocity ±1. Upon collision step pairs annihilate each other. The lower 
lines h-i, h-2, ■ ■ ■ follow the same dynamics, independently of each other, under the condition 
that nucleations are suppressed whenever they violate the monotonicity constraint, i.e. for 
height line i, i < 0, nucleations occur only in the region {x E M| |x| < t, hi^i{x,t) — hi{x,t) > 
2} with space-time rate 2. As for the RSK dynamics, we have to convince ourselves that the 
GW dynamics lives on P^ with probability one. The simple procedure is to use duality. We 
regard the multi-layer PNG model as the stochastic evolution of the random field ?7f(j, x) over 
Z X M with values in {0, 1} by setting 

j 1, if hi{x, t) = j for some i, 
rit[j,x) = < (2.1) 
I 0, otherwise. 

The starting configuration is ?7o(j, a^) = 1 for j < 0, tjqIJjx) = for j > 1. Let us introduce 
the fiipped configuration through f]t{j,x) = 1 — rit{j,x) with the corresponding height lines 
hi{x, t), £ = 1, 2, . . . . Then /ii(x, t) evolves again as the PNG droplet, now growing downwards 
towards negative j. Thus for the original process, at any given time t, there is a random index 
I such that hi{x,t) = i for all x. This proves that at any time t with probability 1 there are 
only finitely many up- and down-steps. 

Proposition 2.2 Let {hi{x,t), x G M, £ G N_} be the height lines as generated by the GW 
dynamics. Under the step map S their joint distribution is fit, fit of Theorem \2. 1\ . 



Proof: We did not discover a global version comparable to the proof of Theorem |2.1| and have 
to rely on an infinitesimal argument. Let L{t) be the forward generator of the Markov jump 
process induced on P^ through the GW dynamics. We have to prove that 

j^{e-'''wt)=me-'''wt (2.2) 
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which means 



-Uwt = L{t)wt , (2.3) 

since dwt/dt = 0. L{t) has four pieces, (i): There is a loss term —At due to nucleation events 
at hQ{x,t) for x G [— which cancels the left hand side of ( |2.3| ). (ii): The free flow terms 
d/dyj'~^ of the up-step motion and —d/dy^^~ of the down-step motion vanish when acting on 
the Lebesgue measure. There is no boundary contribution at ±t, because the boundary moves 
with the same speed as the steps, (iii) and (iv): Let = \{x E R| \hi{x,t) — hi_i{x,t)\ > 2}| 
with I ■ I denoting the one-dimensional Lebesgue measure and let r = Yle<o''^^- -^^^ given |n| 
the current configuration gains in probability due to a transition from |n| + 1 to |n| through 
the collision of an up-step and a down-step. Their relative velocity is 2. Thus for a small time 
interval dt the gain is 2r dt, since the ratio of the weight at |n| + 1 to the weight at |n| equals 1. 
There is a loss of probability due to nucleations in the current configuration. For the time dt 
it is 2r dt, since nucleation events have intensity 2. Thus the gain term (iii) cancels exactly the 
loss term (iv). More extended versions of our argument can be found in |] where, however, 
a different geometry is discussed. □ 

Next we introduce a discrete time version of the multi-layer PNG. As explained in f^, |10 



this model is isomorphic to the shuffling algorithm for the Aztec diamond. We discretize time, 
now denoted by r = 0, 1, . . . . We also discretize the space axis in units of 6. As before, at 
time r, the height lines are h^{x,T), x G M, ^ G N. h^{x,T) = £ for |x| > 6t. The non-crossing 
constraint /i£_i(x, r) < hi{x,T) is in force. Up- and down-steps are allowed only at midpoints 
of the form (m + |)5, m G Z and their distances must be odd, i.e. of the form (2m + 1)5. To 
update from time r to time r + 1 only changes inside the strip [— (r + 1)6, (r + 1)6] are allowed. 
The actual update consists of a deterministic and a stochastic step. 

(i) Deterministic step: Given hi{x,T) every up-step is moved a (5-unit to the left, every down- 
step a (5-unit to the right. If at time r there is a block of length 26, short a 2(5-block, with a 
down-step to the left of an up-step, then they annihilate each other, i.e. in this block hi{x, r) 
is replaced by its maximum. The configuration after the deterministic step is denoted by 
hi{x,T + 1). 

(ii) Stochastic step: The constant pieces of each height line he{x, t + 1) are subdivided in 
consecutive 2(5-blocks. To fix their location, the left endpoint of a 25-block either coincides 
with the right endpoint of a 25-block or is ^6 away from an up-step, resp. from a down-step. 
If ho{x, r + 1) = 0, the 25-blocks are of the form [(-r - 1 + 2m)6, (-r - 1 + 2m + 2)6]. If 
h£{x, T + 1) = i for all x and hij^i{x, r + 1) 7^ £ + 1 for some x, the 25-blocks of hi{x, r + 1) 
start at (?/+ |)5 with y the position of the first up-step (from the left) of hi^i{x, r + 1). Finally 
we disregard those 25-blocks for which h^^i{x, r + 1) — hi{x, r + 1) = 1 for some x inside the 
block. After these preparations the stochastic update can be performed. Independently for 
each 2(5-block, we keep the original piece of the height line with probability 1 — g, 0<g<l, 
and otherwise nucleate an up-step to the left and a down-step to the right midpoint of the two 
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adjacent (5-intervals. The line configuration after tlie stocliastic update is denoted by h^i^x, r+1). 

In tlie limit of rare events the discrete multi-layer PNG model converges to its continuous 
time cousin. We set t = t6 and denote by [t] the integer part of t. Then space-time is discretized 
in cells of lattice spacing 6. A nucleation event covers a block of two adjacent cells. If we set 
q = 4(5^, then in the limit 5 ^ we obtain a Poisson process of intensity 2. Therefore in this 
limit he{x, [t/6]) — >• he{x,t) stochastic process. 

The discretized multi- layer PNG model inherits the simplicity of the time r measure, /i(r). 
The height line he^x, r) has rii up-steps. The total number of up-steps is then X]^<o'^^ ~ 
the collection of height lines {/i^(x,r), a; G M, £ < 0} we assign the weight (g/(l — q))^- The 
partition function, Z^{t), is the sum over all weights. We set Z^{t)~^ = F(^{ho{x,0) = 0}). 
Therefore 

T 

Ur) = 11(1 - q)-^ = (1 - . (2.4) 

If the weight at time r is denoted by w{t), we claim that 

/i(r) = Za{t)-'w{t) (2.5) 

is the time r probability measure of the discrete multi-layer PNG. Let be the transition 
kernel from r to r + 1, as explained in steps (i), (ii) above. We have to show //(r + 1) = /i(r)-R'^, 
equivalently 

{l-qy-^^w{T + l) =w{T)Kr- (2.6) 



is established in Proposition ^]3| below. But first we want to convince ourselves that 
/i(r) yields fit of the continuous time PNG in the limit 6^0. We note that for r = [t/6],q = 
452, 

limZd([t/5]) = lim exp ( - ^[t/5]{[t/5] + 1) log(l -4:5^))= e^*' =Z{t). (2.7) 

A configuration with n up/down-step pairs has the weight (g/1 — g)" = (45^)". Because of 
the constraint in the up-step locations, in the limit 6^0 the weight converges to the 2n- 
dimensional Lebesgue measure constrained to rt(n), n = |n|. Thus fi{[t/5]) fit as 6 0, as 
it should be. 

Proposition 2.3 Let the weight, w{t), of the height lines of the discrete multi-layer PNG be 
given by {q/{l — q))"', where n is the total number of up-steps (equivalently down-steps). Then 
holds. 

Proof. Let w{t + \) be the weight for the configuration /^^(a;, r + 1). We construct from it the 
configuration /i^(x, r + 1) by removing all spikes from /i^(x, r + 1), i.e. all 25-blocks containing an 
up-step to the left and a down-step to the right. Let Si be the number of spikes for /i^+i(x, r + 1). 
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let be the number of up-steps for hi{x,T + 1), and let bg be the number of 25-blocks with 
hi{x, r + 1) constant such that hi^i{x, T+l) — hi{x, r+1) > 2 within that block, Bq is the number 
of flat 2(5-blocks of hQ{x, r + 1). Next we map the configuration hg{x, r + 1) to the configuration 
hg^x, r) by moving all up-down-steps one step backwards in time. By construction, does not 
change. Let ag be the number of downwards open blocks, i.e. flat 2(5-blocks of hi{x^T) such 
that hi_i{x,T) has distance > 2 within that block. 

The transition kernel K^^ = Kf^Kg, where i^d is the deterministic step (i) and Kg is the 
stochastic step (ii). To compute w{T)Kr, we first evaluate w{T)Kd in the configuration he{x, r + 
1). We have to sum over all line configurations hi{x,T) leading to hi{x,T + 1) under K,i. A 
downwards open block of hi{x,T) had either no steps, weight 1, or a downwards spike, weight 
g/(l — q). Summing over these 2"'' possibilities results in the weight (1 — q)~"''^{q/{^ — q))"'' 
for he{x,T + 1). Applying the stochastic transition Kg yields the weight w{T)Kr evaluated at 
hi{x, r + 1) as 

q^^(l-qf^-^^(l-q)-^[q/(l-q)Y^. (2.8) 

On the other hand, according to i/;(r+l), /i^(x,r+l) has the weight {q/ {^—q)Y^^^^ ■ Comparing 
with (|2.6|) and ( p^) , we have to prove 

Y,{h - a,) = r + 1 . (2.9) 

Let be the index of the last height line for which /iAr(x, r) = N for all x. Then oat = 0. 
For two adjacent lines it is easily verified that 

ht- ai+i=ni+i-ni, i < -1 . (2.10) 

Inserting in the left side of (pl9|) and using = 0, the telescoping sum gives 60 + ^^0, which by 
definition is independent of ho{x, r + 1) and equals r + 1. □ 

3 1 + 1-dimensional Fermi field 

Wt is the uniform distribution on all allowed line configurations of the continuous time multi- 
layer PNG. Except for exclusion (entropic repulsion), the height lines do not interact. Such a 
statistical mechanics system is most conveniently analyzed through the transfer matrix method. 
Its implementation requires the height lines to be restricted to a bounded interval {— M, — M + 
1, . . . , M} = Im- The case of interest is then obtained in the limit as M ^ 00. To explain 
the principle, we omit the argument t and label the height lines more conventionally as hi{x), 
i = 1, . . . , N, N < 2M + 1, with x G [0, t]. The height lines are constrained through —M < 
hi{x) < ■ ■ ■ < hi\i{x) < M for all x G [0, t]. In addition we fix the initial configuration q and the 
final configuration q', i.e. /i^(0) = g^, hiit) = q'^, i = 1, . . . , N. As before under the step map 
the i-th height line is specified by the location of the up-steps, < < ■ ■ ■ < y^^ < t, and 
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the location of the down-steps, < ~ < ■ ■ ■ < y'^ < t, where rii 7^ is allowed. Admissible 
line configurations are assumed to have a uniform weight, which means that a small volume 
element has the weight Y[t=i YYjLi Y[^''=i dy^j^ dy^-,~ ■ The configuration with no steps has weight 
1. 

We want to compute the partition function q') which is defined as the weight integrated 
over all admissible step configurations. For this purpose the simplex Vti^ = {q E TL^\ — M < 
gi < ■ ■ ■ < gTV < M} is introduced. Clearly q,q' E VLn and we regard Zt{q,q') as a \^n\ x I^a^I 
matrix. By the product property of the Lebesgue measure Zt satisfies the semigroup property 

ZtZs = Zt+s, t,s >0, Zo = l, (3.1) 

1 the identity matrix. Thus there exists an infinitesimal generator G^, such that 

Zt = e"*^^, t > 0. (3.2) 

Differentiating at t = one concludes that Gj^ acting on functions / on fi^v is given by 

GNf{q) = - Yl c(g,g )/(?'), (3.3) 



where for q,q' G Q 



N 

= ( J = (3.4) 

I otherwise. 

Computationally much more powerful is to impose the constraint of no overlap through 
antisymmetry. Let J-'n be the subspace of ^2((-^m)^) consisting of antisymmetric functions 
over {Im)^ , i-e. / G J^n satisfies 

/(gi, ...,qN) = i-ir'"'' fiq.il), • • • , g.(7V)) (3.5) 
for every permutation n. J^n is equipped with the canonical basis /g, g G VLn, defined through 

* TT 

with 5q{q') = 1 ii q = q' and 5q{q') = otherwise. The normalization is chosen such that 
{fq,fq') = ^q{q'), wlth (■,■) dcuotiug the scalar product in JF^. Let us also define the one- 
particle Hamiltonian 

H^'ij{-M) = -ij{-M + 1), H^'^iM) = -ij{M - 1), 

H^'iP{n) = -'ip{n + l) -ip{n-l) for |n| < M, (3.7) 

as acting on functions tp over Im- The corresponding A^-particle Hamiltonian, is then given 
through 

N 

= 5^ 1 ® ■ ■ ■ ® i/f ® ■ ■ ■ ® 1, (3.8) 
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where H^^ is inserted at the j-th position of the A^-fold product. Clearly for / G one has 
H^^j^f G JFjv and H^^j^ is regarded as acting on jFjy. With these notations one has the identity 

iU e-*<-/,,) = e-*«-(g, q') = Z,{q, q') (3.9) 

for q, q' E ^In- 

At this point it is more convenient to switch to fermionic language which is devised precisely 
to take the antisymmetry into account. The CAR algebra over Im is generated by a*{j), a{j), 
j G Im- They satisfy the canonical anticommutation relations 

{a(z),a*(j)} = 5.,, {a{i),a{j)} = 0, {a*(z), a*(j)} = 0, (3.10) 

i,j = —M, . . . , M, {A, B} = AB + BA. In the Fock representation the algebra is realized as 
operators on the antisymmetric Fock space JF over Im-, 

2M+1 



^= (3.11^ 



N=0 



The second quantization of the (2M + 1) x (2M + 1) matrix is defined by 

H^'= E (3.12) 

H^^ restricted to jF/v agrees with H^^j^. From ( |3.9| ) one concludes that as a fermionic operator 
the transfer matrix is given through 



e-*^", t>o, (3.13) 

which covers all < iV < 2M + 1. 

We exploit the new flexibility by assuming that g, g' G = IJ^Iq^ ^^n, which is identified 
with {0, l}-^'^^ The case of interest is g = g' and each boundary configuration g has the weight 
exp (A(j) In other words the configuration g has a product weight with 

factor e'^^-'^ if site j is occupied and factor 1 if the site j is empty. The corresponding partition 
function is then given through 

Z,\=tr[e^"e-*^d"], (3.14) 

where the trace is over JF and 

N'' =Y,Hj)a*ij)a{j). (3.15) 

The probability that site j at x is occupied, i.e. hi{x) = j for some i, is obtained from the 
transfer matrix as 

(Z,\)^^tr[e^^^e-^'^^^^a*(j)a(j)e-(*-^-)^^"] = tr[pf a*(j)a(j)], (3.16) 
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with the density matrix 



P: 



-M 



(3.17) 



< X < t. 

Exponentials of operators quadratic in a, a* are easily handled. If A is a (2M + 1) x [2M + 1] 
matrix with second quantization 



A= ^ a*{i)Aija{j), 



(3.18) 



then 



on jF/v. This implies 



tr[e^] = det(l + e^) = Zj 



(3.19) 



(3.20) 



compare with [|TI|. Let us set pa = (^a) as density matrix. The two-point function has 
the form 



R{t,j) = tr[pAa*{t)a{j)] = (^(1 + e 
and more generally 

tY[pAa*{ii) ■ ■ ■a*{im)a{in) ■ ■ ■a(ji)] = ^m,n det {R{ik, jk')) 



(3.21) 



(3.22) 



The expectations of other monomials are determined by means of the anticommutation relations 
( p.lOp . One may regard tr[p^-] = uja{-) as a linear functional on the CAR algebra. By 



definition uja{^) = 1- If in addition uja is positive, then uia is called a quasifree state fTS 



As can be seen from (|3.19|) products of exponentials follow the same pattern. For the 
(2M + 1) X (2M + 1) matrices A, B we set 



Then 



(3.23) 



(3.24) 



with ^ defined as in ( p. 181) . 

Each a*{i)a{i) is a symmetric projection and thus has eigenvalues in {0, 1}. {a*(z)a(z), i G 
Im} is a family of commuting operators. Under p*^ they have a joint signed spectral measure 
which by construction is a probability measure on {0,1}^^^ By (|]22D it is of determinantal 
form, compare also with ( |3.54| ) below. Thus under the family {a* (2)0(2), i G Im} is a point 
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process on Jjv/. In the probabilistic literature point processes of this structure are known as 
determinantal |T^]. 



Our construction extends to occupation variables depending on several x. For example, the 
probability that site j at x and site iaXy,0<x<y<t, are both occupied is obtained through 
the transfer matrix as 

(Z^,)-'tr[e^''e-^^d''a*(j>(j>-(^-^)^"a*(z)a(z)e-(*-^)^d'']^ (3.25) 

Such expressions can be computed using ( |3.24|) and (|3.23| ), and the rules 



a{j)e-'^ = E e*V(j)e-*^ = «*(0(e*^).r (3-26) 



With these preparations we return to the PNG droplet. Recall the definition ( p.l|) of the 
random field of occupation variables rit{j,x), j G Z, x G M: rit{j,x) = 1 if he{x,t) = j for some 
i and T]t{j,x) = otherwise. By definition T]t{j,x) = for j > 1, |x| > t, and r]t{j,x) = 1 
for j < 0, |x| > t. The joint distribution of T]t{j,x), j G Z, x G M, is induced through the 
probability measure fit with expectation denoted by Ej. Our first goal is to obtain the joint 
distribution of f]tij,0), j G Z. From the considerations above, it is obvious that this point 
process is determinantal. The only remaining task is to compute the two-point function and to 
study its limit behavior. One should pay attention to a minor linguistic problem. The meaning 
of the time parameter t of the PNG model is now reduced to a mere scaling parameter, j labels 
fermionic space and x stands for fermionic time. Thus space-time means now Z x M. 

Let us start by explaining the result for moments at x = 0. We introduce the limit M — > oo 
of the one-particle Hamiltonian in ( p.7|) as 



H^tP^n) = -ij{n + l) -ij{n-l). (3.27) 
In addition, we add a linear potential of slope 1/t to define 

n 

Htip{n) = -ip{n + 1) - V^(n - 1) + --^{n), (3.28) 

regarded as an operator on £2 = •^2(Z). Ht has the complete set of eigenfunctions Lp^''\n) = 
Jn~.i{2t), I G Z, with eigenvalues Si = ^, Ht^p'^^'' = eiip^''\ {ip'^^\ ip'^'-">) = dw, (-, ■) denoting the 
scalar product. Here Jn{z) is the Bessel function of integer order n and we follow throughout 
the conventions of [jl4[. Chapter 9. We will need the spectral projection Bt onto {Ht < 0}. In 
position space its integral kernel is the discrete Bessel kernel 

Btit.j) = J'W^t)J,.i{2t). (3.29) 
z<o 



Using that Ht^p^''^ = £i^^^\ (|3.29|) can be converted into a telescoping sum with the result 



Btii,j) = -^(J.-i(2t)J,(2t) - J,(2t)J,_i(2t)) (3.30) 
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for i ^ j and on the diagonal 



Bt{i,i) = t{Li^i{2t)Ji{2t) - Li(2t)J,_i(2t)), (3.31) 

where Lj{2t) = ^Jj(2t). We also introduce the CAR algebra over Z. It is generated by the 
operators a{j),a*{j), j G Z, satisfying the canonical anticommutation relations (|3.1CI|) . Let Ut 
be the quasifree state on the CAR algebra defined through u;t(a(j)) = = LJt[a*{j)^ and 
the two-point function 

u;t{a*{i)a{j))=BS,j), (3.32) 

which means that higher order monomials satisfy (|3.22| ) with R replaced by Bt [|l^. Ut is the 
ground state for non-interacting fermions with one-particle Hamiltonian ( |3.28D at zero chemical 
potential. 

Theorem 3.1 We have 

Et(f['nt{jk,0)] = u;t(f[a*{jk)a{jk)] = det {Bt{jk,Jk')),^,^,,^^, (3.33) 
the second equality being valid only for pairwise distinct points ji, . . . ,jm- 

Proof: The state Ut has an infinite number of fermions and cannot be represented as a vector in 
Fock space. Thus we first have to constrain to finite volume Im, such that \he{x,t)\ < M,i = 
0, . . . , —M. The corresponding uniform distribution is denoted by /if^. Clearly fif^ converges 
to fit as M oo. By construction hi(t,t) = hi{—t,t) = i, £ = 0, —1, —2, ... . In other words 
Vt{j^ ~t) = Vtijjt) = 1 for j < and = for j > 1. It is convenient to approximate these 
boundary configurations through the weight e^, — M < j < 0, weight e~^, 1 < j < M, for an 
occupied site and weight 1 for an empty site in the limit (3 —>■ oo. Thus, if we set A^*^ as in 
(|1|) with A(j) = 1 for -M < j < and A(j) = -1 for 1 < j < M, we have 

n r/*(j., 0) ) = lim ^tr[e^^^^e-^.^^ H ^^U^Hj'^y-'''^] (3-34) 
fc=i ^ k=i 

with Z{(3) the normalizing partition function. Since the moments of a quasifree state are 
determined by the two-point function, to prove (|3.33| ) it suffices to consider the expectation 
of a*{j)a{i) and to subsequently take the limit M oo. Let P^^ be the projection onto 
{-M, . . . , 0}, P^^ onto {1, . . . , M}, Pj^'^ + P*'^ being the identity. Then 

hm Z(/5)-Hr[e'3^"e-*^"a*(j>(z)e-*^"] 
= hrn^ ((1 + e*^"(e^Pf + e-^P^^)e*^")-i 

= (e-*^"P^-^(Pf + pM^-2tH^'pMylpM^-tHi^^ (3 35) 
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If P_ denotes the projection onto N_, P+ = 1 — P_, then 

hm = P±, hm e'"^' = e*^^ (3.36) 
To prove the theorem we only have to check the identity 

e-*-^^P_(P+ + P_e~^*^^'P_)"^P_e-*^<^ = Bt (3.37) 

as an operator identity on £2- 

We define the left shift D, Dipiji) = ■il){n + 1), and the adjoint right shift D* , D*il){n) = 
^{n - 1). Clearly [D,D*] = 0. One has Hd = -D- D*. Using 2|J„(t) = J„_i(t) - Jn+i{t), 
one obtains 

j^Bt = (D* - D)Bt - BtiD* - D) . (3.38) 
Integrating with the initial condition Bq = P_ yields 

and 

Therefore ( p.37| ) is equivalent to 

^-2tDp_^-2tD* ^ ^p_^2tD*^2tDp_yl ^3 ^^) 

as an operator identity on P_£2- We decompose our space as £2 = P+C-2 ® P--^2- Then with the 
definition 

^-2tD =h ^\ (3.42) 



h c 



we have 



Using the splitting of (|3.42|), one constructs the inverse operators e^*^,e^*^ . By a straightfor- 



ward computation one obtains 



(c*c) 

From Theorem one immediately infers the distribution of the height of the PNG droplet 
at X = 0. Clearly 

Pi({/i(0,t) < n}) = Pi ({r/j(j, 0) =0 for all J > n}) 

= lim uoti fTe-^'^*(^>(^)^ = lim det (l - (1 - e-P)PnBt) 

^ j=n 

= det(l-P„Pt), (3.45) 
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where P„ denotes the projection onto {n, n + 1, . . . } in £2- Since L(0, t) = h{0, t), see eq. (|1.8|) , 
we have rederived that the length of the longest increasing subsequence of a Poissonized random 
permutation has a distribution linked to the discrete Bessel kernel. Previous proofs take the 
route via the Plancherel measure. We refer to [|l6l. It would be of interest to better understand 
how these proofs are linked to the multi-layer PNG. 

So far we considered only the distribution of rit{j,x) at x = 0. The transfer matrix method 
can handle also the distribution referring to several x, like the joint distribution of rjtihO), 
rit{j,x), see eq. ( |3.25| ). The transfer matrix is generated by the Hamiltonian of ( |3.27| ). As 
in the case of fixed x, the joint moments have determinantal form with the entries given by the 
Euclidean Fermi propagator Bt{j,x]j',x'). Following the scheme in (|3.34|) it is defined through 
a finite volume approximation, Bt{j,x] j' ,x') = linijv/^oo ^^"^/s^oo Bf'^'^ljjX; j' ,x'), where 

for —t<x<x'<t, 

for —t<x'<x<t. 

(3.46) 

Note that time order must be respected in such a way that there are only decaying exponentials. 
The minus sign in ( |3.46| ) for x' < x results from commuting a* and a. At coinciding arguments 
the definition conforms with E(?7t(j, x)) = Bt{j,x;j,x). Using ( p.26|) one obtains 



Bf^U^x^i'.x) = < 



B,ij,x;f,x') = (e"^^^(i?, - M(x - x'))e^-'^^) . ., 



(3.47) 



for |x| < t, \x'\ < t, X x', with the step function 6{x) = for x < 0, 6{x) = 1 for x > 0. B^ 
has a jump discontinuity at x = x'. For coinciding time arguments one has 



Bt{j,x;j',x) 



IJ3 



(3.48) 



For later use the propagator is rewritten in the eigenbasis of Ht. The integer order Bessel 
function has the representation 



1 



Jn{2t) = — —e 
Ztii 



where the contour integration is a circle around 2; = 0. Therefore 



dz 



(3.49) 



(3.50) 



Substituting z by {t + x)^/^(t — xY^'^z yields 



(e-^^V.(2t))^ = J42v/t^^) ( 



t + X\'^/2 



t — X 



(3.51) 
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and thus for x ^ x' 

= Y,^gn{x' -x)9{{x-x'){l + \))[j^) J,.,{2VW^ 



xJ,,_,(2v/^^)(^)"'"''^ (3.52) 



At coinciding arguments x = x' one has 

't + X\(i-j')/2 



B,{j,x;j',x) = J2 (^) J,-i{2VW^J,,.i{2VW^ (3.53) 

KO 

With these preparations for a general moment of the density field rit{j,x) one has the identity 

(m 
k=l 

= hm hm Z{(3)-hi\e^^" e-'^"^' ( IT e-^-(^)^"a*(j^(M)a(j^(fc))e^--w^d' V"*^"" 

M— >oo /3— >oo L ^fc ' ' 

= det (Bt{jk,Xk]jk',Xk')) . (3.54) 

V / l<k,k'<m 



As written, ( p. 54]) is valid only for pairwise distinct Xi, . . . , Xm, where tt is the unique permu- 
tation of 1, . . . ,m such that the time ordering —t < Xn{i) < ■ ■ ■ < < t is ensured. The 
spatial arguments ji, ■ . . ,jm are arbitrary. While each off-diagonal factor in the determinant 
has a jump discontinuity at Xk = xy-, the determinant itself depends continuously on xi, . . . , x^, 
and thereby the continuous extension of ( p.54| ) holds for all G Z, — t < < t, = 1, . . . , m. 



As an application of (|3.54|) we establish the joint distribution of {rit{j,x),j G Z} for fixed 



X, |x| < t. From ( |3.53|) one derives immediately 



where g is a multiplication operator with diagonal entries g{n) = ((t + x)/(t — x))"^^. In (|3.54| ) 
we take the limit of coinciding Xk = x, k = 1, . . . ,m, leaving ji, . . . , pairwise distinct. Upon 
forming the determinant in ( p.54D the similarity transformation g drops out and the result is 



( p.33| ) with Bt replaced by Bnp—^. Thus the joint distribution of {rit{j,x),j G Z} is again 



given through the discrete Bessel kernel with time parameter modified from t to \^t'^ — x^. 

The same conclusion can be drawn by taking the discrete PNG model as starting point. 
As explained in the analogue of the fixed x distributions is given through the Krawtchouk 
polynomials. Their limit as 5 — > 0, g = 45^, 6t = t, yields the joint distribution of {f]tij, x),j G 
Z} as given through ( p.35 ) with parameter \Jt'^ — x^ instead of t. 



In the following section we establish the scaling limit of the PNG droplet at locations of 
order {%)t^/^,2t + ut^/^), y,u eM.. Since at x = wt, \w\ < 1, the distribution is determined by 
the discrete Bessel kernel B^^ ^_^i , we could instead of = choose any other reference point 
(wt, 2\J\ — wH) and relative displacements {wt + yt^^^, 2\J\ — wH + nt^/'^). Except for scale 
factors, the limit t ^ oo does not depend on the choice of w. 



18 



4 Edge scaling, convergence to the Airy process 



We plan to establish that the statistics of the PNG droplet close to x = converges to the Airy 
process. Since only moments are under control, the natural strategy is to prove that the Fermi 
field of Section § has a limit when viewed from the density edge. In particular, this implies, 
that the statistics of the last fermionic world line has a limit, which is the desired result. 

Since {h{0,t)) = 2t, the focus has to be at x = + yt"', j = 2t + ut^, y,u & M. fixed. 
Rescaling dip/dx = Hfip accordingly one obtains 

= r ( - ^(m + t-^) - tlj{u - r^) + ^(2t + ut^)ij{u)) . (4.1) 

To have a limit operator as t — > oo requires 

« = | 13 = I, (4.2) 
and with this choice ( [4.1D converges to 

regarded as a self-adjoint operator on ^^(IR). H is the Airy operator. The limit density field 
must correspond to free fermions with H as one-particle Hamiltonian. The fermions are in their 
ground state at zero chemical potential. 

Let us first describe the Fermi field in more detail. The Airy operator H has M as spectrum, 
which is purely absolutely continuous. The generalized eigenfunctions are the Airy functions, 

Ai(n - A) + nAi(n - A) = AAi(n - A) . (4.4) 



In particular the completeness relation 



dX Ai{u - X)Ai{v - A) = 5{u - v) (4.5) 
holds. K denotes the spectral projection onto {H < 0}. Its integral kernel is the Airy kernel 



K{u,v) = I dXAi{u- X)Ai{v - X) 

-oo 

1 



U — V 



(Ai(n)Ai'(iO - Ai'{u)Ai{v)) . (4.6) 



Next we introduce the Fermi field a{u), a*{u), indexed by n G M. To distinguish from the 
fermions on a lattice we should use a different symbol. Since the latter will not reappear, we 
find it more convenient to stick to familiar notation. Integrated over a test function / G L2(M) 
the Fermi field becomes a(/) = / du f*{u)a{u), a*(/) = / du f{u)a*{u) = a{f)*. They satisfy 
the canonical anticommutation relations {ci{f) , a{g)} = = {a*{f),a*{g)} and {a{f),a*{g)} = 
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{f,g) with (-, ■) denoting the inner product of / G L^(M) |T^, and generate the CAR algebra 
A. On A we define the quasifree state u through a;(a(/)) = = Ci;(a*(/)) and uj(^a*{f)a{g)^ = 
(/, Kg). In particular, the moments of the density field are given by 



uji Y\a*{uk)a{uk)j = det K{uk,Uk')i<k,k'<m, (4.7) 



for pairwise distinct mi, . . . compare with (|3.33|) . ( [4.7|) is the m-th correlation function. It 
vanishes at coinciding points. 

To extend to unequal times one defines the Euclidean propagator 

K{u,y;u\y') = (e"^^ (A' - Iie(y - y'))e^'^) (n, n') 

sign(?/' -y) [ d\e{\{y - y'))e^^y'~y^Ai{u - \)Ai{v - A), (4.8) 



for y ^ y\ written in terms of eigenfunctions of iJ, and 

K{u,y-u,y) = K{u,u), (4.9) 



compare with ( p.4?| ), (|3.48|) . The quasifree state uj and the propagator are both determined by 



the Airy operator if, which implies that K depends only on y — y'. 

The Airy field, denoted by ^(/, y) = f du f{u)^{u, y), is the density field of the Fermi system 
defined through ( [4.8|) , ([4.9|) . As in (|3.54|) , its moments are of determinantal form and given by 

k=l ' k=l 



As it stands the left hand side of ([4.10|) is only defined for pairwise distinct yi, . . . ,ym, but as 
in (|3.54|) it can be continuously extended to arbitrary time arguments. Since K depends only 



on y — y', the Airy field ^(/, y) is stationary in y. 

Having introduced the limit object we turn to the edge scaling. Recall that the PNG droplet 
has curvature. Therefore we set the scaled density field of the multi-layer PNG model, denoted 
by ^j, as 

Uu,y) =t'/\{[2t + t'/'{u-y%t'/'y) , (4.11) 

[■] denoting the integer part. When integrated over the real, smooth, and rapidly decreasing 
test function / we have 



^t{f,y) = J duf{u)it{u,y) 

oo 

= E fit~'/'{j - 2t) + y')vt{j,t'/^) + 0(t-V3) 

j=-oo 

= 5^/(r^/'(Mt'/'2/,^)-2t) + 2/')+0(ri/3), (4.12) 
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where the error of order t results from integrating over cells of size t in the defining 
identity. Thus (4.12) shows that through controlling the limiting moments of ^t{I\y) one can 



infer the limit of the scaled height lines t ^/^{hi{iP'/^y,t) — 2t) + y^. 

Theorem 4.1 Let fi,-.-,fm be smooth test functions of compact support. Then the following 
limit holds, 

\imEJf[Ufk,yk)) = E(f[af,,y,)] . (4.13) 



Proof: Comparing ( ^.54[ ) and ( [4. 101 ) we have to establish that the propagator ( |3.47| ), properly 



scaled, converges to the continuum propagator ( |4.8| ). This limit can be handled most directly 
in the representation ( p.52 ). We will need a separate argument for y ^ y' and for the left, resp. 



right, limit y = y'. 

The case y < y' runs in complete parallel to y > y'. To simplify notation let us assume 
y < y' . The propagator for the scaled density field is 

Uu, y- u\ y') = e-'''''ye^^~y'H'/'B,{[2t + e'\u - y%t'''y, [2t + t^l\u' - y'%t'/'y') 

y, y' are fixed and Kt is considered as a function on M^. It is constant over cells of size t^^/^. 
We used here the freedom that the determinant of (|3.54|) does not change under a similarity 
transformation and multiplied with the factor exp (2t'^^^{y' — y)) exp ((m — y'^)y — [u' — y'^)y') 
which diverges a.s t —>■ oo. Kt, with the obvious extension to y > y', determines the moments 
of the scaled density field as 

(rrt \ „ m 

l[Ufk,yk)] = / l[dukfk{uk) det {Kt{uk,yk;uk',yk')),^,,,^^ + 0{t-'/^), (4.15) 
k=l ^ k=l 



clearly analoguous to ( [4.10[) 



We insert ( P32|) into Qm] ). Then 



Kt{u,y;u',y') = t'^'' Yl e^''-'^'t'/' JlMfHu-y^-i)] {^Wl - t'^'y') 



i6i-i/3r 



{2t^l-t^/y^) 



X <; exp {-{2t^/'y - (u - y')y - ly)) { j^^/g 



l+t-'i'y 

y 



1/3 N (2t+tl/3(„_y2_;))/2 



X exp ((2t^/V - (- - y V - V)) 



(4.16) 
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\i y < y' and / < 0, the term { ■ ■ ■ } is uniformly bounded in t, / and converges to 1 as t — oo. 
In fact this holds uniformly in u, u' on compact sets. By a result of Landau [|T^ 



sup 



(4.17) 



with c = 0.7857- • •. From the asymptotics of integer Bessel functions, cf. (9.3.23) of [14 
conclude, uniformly for u varying over a compact set. 



we 



^hm Jy2t+tU^iu-y^)] (2tV 1 - t2/3y2) = Ai( 

Since e^^'~^^', / < 0, is integrable, by dominated convergence 

lim Ktiu, y; u', y') = K{u, y; u', y') 



u 



f— >oo 



(4.18) 



(4.19) 



uniformly over compact u, u' sets. 

Next we consider y' \ y in (|4.14]) , the right hand limit being handled analogously. If the 
discrete Bessel kernel is transformed according to ( 3.55|) , then 



Kt{u,y;v,y) 



1 - t~y^y 



1 + t-i/3y 



X 



^/^3p^([2t + t'/\u - y% [2t + t"\v - y')]) . 



(4.20) 



The first factor is uniformly bounded over compact u,v sets and converges to 1 as t — oo. By 
Proposition (4.1) of ||18[ the discrete Bessel kernel with our scaling converges to the Airy kernel 
K(u, v) uniformly on compact m, v sets as t ^ oo. 

We conclude that ( [4.19| ) holds not only for y ^ y' but also for its right and left limits. Our 
claim follows by taking the limit t — > oo in ([4.15|) which then yields ( [4.10| ). □ 

The Airy field .^(/, y) is stationary in y. ^{f, y) is a point process for fixed y. Its average 
density is given by 



E(^(m, y)) = -u Ai{uY + Ai'(M) 



(4.21) 



which has the asymptotics 

E(e(w,y)) - 



i^cos(4|M|3/V3) + C(|M|-5/2) for 



47r|M| 



u 



-oo, 



967r 



exp ( - 4^3/2/3) 



for M — >• oo. 



(4.22) 



Note that for m — oo the density decays quickly because of the increasing linear potential, 
whereas for n ^ — oo the density is limited through the Fermi exclusion. In particular, the 
point process for rj{f\y) has a last point at hQ^y) with probability one. Since all points are 
distinct [0, one can label as 



ho{y) > h^i{y) > ... . 



(4.23) 
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y I— hi{y), i G N_, are the fermionic world lines underlying the Airy field. As to be shown 
in Appendix A y hi{y) is continuous with probability one. Moreover {{h^{y) — h£{y')Y) — 
2\y — y'\, which suggests that the path measure for {hi{y), \y\ < c, i = — M, . . . , 0} is absolutely 
continuous with respect to the Wiener measure, i.e. locally hi^y) is a modified Brownian motion. 
Our main focus is the last fermion line. 

Definition 4.2 Let C,{f,y) be the Airy field. The last world line, h^iy), is called the Airy 
process and denoted by A{y). 

We collect the basic properties of the Airy process. 

Theorem 4.3 The Airy process A{y) has continuous sample paths. A{y) is stationary. For 
given y, A{y) has the distribution of X2 of Tracy-Widom, see below eq. ( |7j|j . 

The convergence of the multi-layer PNG model to the Airy field implies that the shape 
fiuctuations of the PNG droplet converge to the Airy process as t — oo. The following theorem 
is the precise version of the main result, Theorem |1 . 1| , stated in the Introduction. 



Theorem 4.4 Let h{x, t) be the height of the PNG droplet and ht{y) its scaled version according 
to ^1.3[ ). Let A{y) be the Airy process. Then for any m, yj, G M, j = 1, . . . , m, we have 

lim Ft{{ht{yj) + < a„ j = 1, . . . , m}) = P({A(y,) < a„ j = 1, . . . , m}) . (4.24) 

Proof. Let fj be the indicator of the interval (oj, oo). Then ( |4.24| ) means 

m m 

^hmP,(f|{6(/„%)=0}) =p(f|{e(/„%) = 0}). (4.25) 

We choose a sufficiently large and split as fj = + g^, where fj is the indicator function 
of the interval {aj,a] and g"" is the one of (a, oo). Then ^t{fj,yj) = (,t{fj',yj) + (,t{g"',yj)- By 
Theorem [4.1| the joint moments of ^t{fj, Uj), j = 1, ■ ■ ■ ,m, converge to their limit. Since their 
limit measure on is uniquely defined by its moments we conclude that ( [4.25|) holds with fj 
replaced by fj. Up to constants the error term is bounded by a sum over terms of the form 

P^(6(^7^%•)>1)<E^(6(^?^%•))= Yl B^^yjzFm^U.j) (4.26) 

which has a bound C{a) uniform in t such that C{a) — exponentially as a — oo, compare 
with (1021 ). □ 



5 Some properties of the Airy process, two-point func- 
tion 

The scale invariant statistics of the PNG droplet is governed by the Airy process. To gain 
some more quantitative information we have to study the Airy process, most prominently its 
distribution at a single point and at two points. 
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We denote by Pa the projection onto the interval (a, oo), i.e. Paipiu) = Xaiu)ip{u) with 
Xa{u) the indicator function of the set (a, oo). By definition 

^{{A{y) < a}) = ¥{{^{xa, y) = 0}) . (5.1) 
Let us set N(xa) = a*{u)a{u)du. Then 



m{Xa,y) = 0) = hm ^(e-^^(^'')). (5.2) 

/3 — ^oo 



Since u is quasifree, 



a; e 



-mixa) 



) =det[l + (e-^-l)Pai^], (5.3) 



where the determinant is in L^(]R). PaK is of trace class [Q and taking j3 ^ oo yields 

F{{A{y) < a}) = det[l - PaK] . (5.4) 

This determinant is studied in and shown to be related to the Painleve II differential equa- 
tion. A plot for the probability distribution of A{y) can be found, for example, in P, ^. 

The next quantity of interest is the joint distribution of A{0), A{y), where by reversibility 
it suffices to consider y > 0. By the same scheme as before one computes 

F{{A{0)<a,A{y)<b}) = P({e(xa, 0) = 0} n {^(x;,, 2/) = 0}) 

= lima.(e-^^(^")e-^^e-^^to)e^^). (5.5) 

Since oj is quasifree, 

^(e-/3iv(x.)g-.Hg-/37V(x,)gyi?) = det(l - 5^) (5.6) 

with 

Bfi = {\- e"^) {PaK + e-^^Pf,e^^i^) - (1 - e"^) V,e-^^P;,e^^i^. (5.7) 
e~y^ P^e^^ K is trace class, cf. Appendix B. Thus 

P({A(0) < a, A{y) < 6}) = det[l - B\ (5.8) 

and 

B = PaK + e-y^P^e^^K - PaC'^^ PbC^^ K. (5.9) 

Clearly, the determinant converges to 1 as a, 6 ^ oo and to as a, 6 — — oo. 

The properties of (|5.4| ) suggest that also the joint distribution might satisfy a differential 
equation. We did not succeed in finding one. Since the main interest is large y, we rely on 
standard asymptotics by employing the expansion 

oo _ 

logdet[l -B] = -J2 -tr[S"] . (5.10) 

n=l ^ 
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When taking the trace of B^, we see that ey^K is a bounded operator, but e ^ remains 
unbalanced, since H is not bounded from below. E.g., tr[i?] diverges as exp(y^/^) for y oo. 
The form tr[i?"] is not suited for studying large y. 

Such a situation is well known in the theory of Fermi systems pO|, |2l||. Since the Dirac sea is 
filled up to energy zero, one has to work with a new representation of the CAR algebra, which 
means to introduce field operators for the particles (energy > 0) and for the holes (energy < 0). 
In this representation the Hamiltonian is positive. There is no need here to enter into the full 



theory. It suffices to note that the series in ( [5. 101) can be resummed such that only decaying 



exponentials appear. Let ct„ = (o"(l), . . . , cr(?T,)) be an n-letter word where each letter is either 
a or b, (Jq = (Jn- Then 



oo ^ n 



logP({A(0) < a,A{y) < b}) = - ^ tr [ J] P.(,)i^.(,v(,+i)] , (5.11) 



"=i W„} j=i 

a{n + 1) = cr(l), with the following convention 



Kh,a = ey^'K. (5.12) 



The large y asymptotics is extracted from ( |5.11 



The spectral representation of e~y^{K — 1), resp. e^^K, yields an integral of the form 
— dXe~'^yg^{X), resp. J^^dXe^yg^{X), with some spectral functions g^, g^. For large 
y the weight concentrates at A = and results in the asymptotics —g{0)/y, resp. g{0)/y. 
Therefore a summand in (|5.11| ) decays as y^", where a is the number of broken bonds, i.e. 



. . .ab . . . and ... 6a ... , in the word (T„ of length n. The only words with no broken bonds are 
either all a's or all 6's. They sum to 

det(l - PaK) det(l - PbK). (5.13) 

As to be expected, the Airy process is mixing and far apart events become independent. 
To order y^"^ we only allow two broken bonds which for large y leads to 

oo oo 



y~' E E ^— (^ai^)'"PaAi)(Ai, (ni^)"nAi) 

m=l n=l 

POO 

/ dX (Ai, P,K(e^ - PaKy'PaAi){Ai, P,K{e^ - P,KY^P,ki) , (5.14) 



Here (Ai, ■ Ai) means inner product in position space with respect to the Airy function Ai(n), 
i.e. for some operator R with integral kernel R{u,v), (Ai, i?Ai) = j du J dvAi{u)R{u,v)Ai{v). 
To compute the two-point function to leading order in 1/y we integrate the probability measure 
with distribution function ( |5.8| ) against a and b and insert to leading order from ( |5.13| ), ( |5.14| ). 
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Using that det(l — PaK) has a good decay for a — oo and (Ai, PaK{e^ — PaK) ^P^Ai) has 
a good decay for a ^ oo, one is allowed to integrate by parts to obtain 



{A{Q)A{y))-{Am{A{y)) = 

dX Jdadet{l-P,K){Ai,PaK{e^-P,K)-^PaAi) + 0(r 



2/^ Jo 



(5.15) 



The Airy process is positively correlated and has a slow decay as l/y"^. 



6 Conclusions 

The height statistics of the PNG droplet, for large t and under the scaling (|1.3|) , are given by 
the Airy process. By universality such a result should be valid for any one-dimensional growth 
model in the KPZ class. The only condition is that locally the macroscopic shape must have 
a non-zero curvature. If the interface is flat on the average, other universal distributions will 
show up 0. In such a situation, at present, no information on the multi-point statistics is 
available. For example in the PNG droplet we could lift the restriction that nucleation events 
are allowed only above the ground layer [— t, t]. By translation invariance of the dynamics 
and the initial condition, h{x,0) = 0, we have {h{x,t)) = 2t for large t and the process 
y t^^/^{h{yt^^^^t) — 2t) = h^{y) is stationary. For fixed y, h^{y) converges to the GOE 
Tracy-Widom distribution P, The problem of the joint distribution of h^{yi), h^{y2) 

remains open. 

The Airy process contains a wealth of statistical information, which cannot be resolved 
easily in Monte-Carlo simulations. In the standard numerical experiment one merely considers 
the second moment of the height differences. It is convenient to subtract the asymptotic mean 
as h{x^ t) = h{x, t) — 2\Jt'^ — x^. The quantity of interest is then 

((/i(x,t)-/i(0,t))') = Gt(x) (6.1) 

for large t. Our main result says that Gt{x) is of scaling form and given by 

Gt{x) ^ t^^^g{t-^/^x). (6.2) 

The scaling function g can be expressed through the two-point function of the Airy process as 

9{y) = {{A{y)-A{0)Y). (6.3) 

For small y, one has 

g{y)^2\y\, y^O. (6.4) 
On the other hand, for large y, A{y) becomes independent from A{0) and 

g{y) ~ 2a2, \y\ oo (6.5) 
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with 02 the truncated second moment of X2; ^2 = (X2) ^ (X2)^ — 0.81320 numerically. The 
asymptotics ( |5.15| ) says that 

- 2a2 - c|?/|"^ \y\^oo, (6.6) 



to next order. Because of the inverse operator in (|5.15| ) the positive constant c cannot be readily 



evaluated. A differential equation, like the Painleve II for the single point distribution, would 
be helpful. 



Appendix A: Continuous sample paths of the Airy process 



Let t I— > X{t) be a stochastic process with values in M. By a criterion of Kolmogorov, cf. [p^] , 
Theorem 2.23, if for some constant c 

E{\X{t)-X{s)\^) <c\t-s\\ (A.l) 

then t (-^ X(t) is continuous (in fact Holder continuous with exponent < |) with probability 
one. Since ( |5.8D provides the joint distribution of A{0), A{y), (|A.1| ) should be an easy exercise. 
We did not succeed and rely on a more indirect argument. 

Let the test function / be smooth and of compact support. By the results of Section ^we 
have 

In the Lemma below we will prove that 

Hm,y)-My')r)<Cfiy-yr- (A.3) 

Therefore y 1— > C,{f,y) is continuous with probability one. Since the vague topology on locally 
finite point measures is countably generated, the trajectory y 1— > ^^<o ^i^^iy) ~'^) continuous 
in the vague topology with probability one. The convergence of a sequence of locally finite point 
measures in the vague topology is equivalent to the convergence of each atom Thus ( |A.3| ) 



implies that y t-^ he{y) for each £ is continuous with probability one. In particular A{y) = ho{y) 
is continuous. 

Lemma A.l Let f be smooth and of compact support. Then there is a constant Cf such that 

E{m,y)-af,0)Y)<Cfy'. (A.4) 

Proof: As a warm-up, and to fix notation, we first compute the second moment. We suppress 
/ and set C,{f, y) = We take y > 0. y < follows from reversal symmetry y —y and 
stationarity in y. As shorthand we define L = K — 1. —L is the projection operator onto H > 0. 
We regard / as a multiplication operator, fip{u) = f{u)tlj{u), and formally set fy = e^^/e~^^. 
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By construction, in the expressions below, only the bounded operators Ke^^ and e L appear 
as factors. By the determinantal formula (4.10) 



= -2{tT{fKfL) - tiifKfyL)) = -2tr(/Al/ - /,)L). (A.5) 

From the spectral representation, one concludes differentiability of any order in y, y > 0, and 
up to linear order 

m^y-^of) = 2ytTifK[HJ]L) + Oiy') 

= 2y{tT{fK[H, f]K) - tr{fK[H, /])) + Oiy'). (A.6) 

For real operators one has tr(AB) = (tr(A5))* = tr(BM*). Since [H, f]* = -[H, /], the first 
summand vanishes and 

m^y-^o?) = ytT{fK[HJ]L) + 0{y')=ytT{K[f,[HJ]]) + 0{y') 

= 2yivK{ff + 0{y^). (A.7) 



The variance ([A.7|) implies that E((/i£(?/) — /^^(O))^) = 2\y\ + for each d in the limit ?/ — ^ 0, 



as to be expected from the construction of the Airy field. 

The fourth moment requires more effort. We have, using stationarity and reversibility in y, 

Ki^y - iof) = 2{mt) - 4E(&) + milil))- (A.8) 



We again use the determinant formula ( [4.10|) , which most conveniently is decomposed into 



cycles. For the fourth moment there are 4! permutations. They subdivide into (i) four 1-cycles 
(1 term), (ii) two 1-cycles plus one 2-cycle (6 terms), (iii) two 2-cycles (3 terms), (iv) one 1-cycle 
plus one 3-cycle (8 terms), and (v) one 4-cycle (6 terms). The sign of the permutation will be 
of no significance for the argument. 

(i) and (ii) vanish, (iii) is the "Gaussian" term, 

(iii) = 3E((e,-eo)')'<c/y2 (A.9) 

from ( |A.7] ). For (iv) we have a 1-cycle and two 3- cycles in reverse order. Summing over all such 
cycles yields 

(iv) = 12tT{Kf){-tT{fKfyLfL)+tT{fKfyLfyL) 

-tiifKfKfyL) + tlifKfyKfyL)). (A.IO) 

Since tiifKfyLfL)* = tifKfyLfyL and tiifKfKfyL)* = tifKfyKfyL the contribution (iv) 
vanishes. 
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It remains to study the 4-cycles. For their sum one obtains 

(v) = tT{fKfKfK{f - fy)L) - Stiif KfK{f - fy)KfL) 
+ti{fKfK{f - fy)LfL) - 3ti{fKfKfLif - fy)L) 
+ii{fKfK{f - fy)LfL) - 3tr(/i^(/ - fy)KfyLfL) 
+tT{fKfLfK{f - fy)L) - 3tT{fK{f - fy)LfKfyL) 
+tr(/i^(/ - fy)LfKfL) - 3tiifKfyLfK{f - fy)L) 
+tr(/K(/ - fy)LfLfL) - 3tT{fKfyL{f - fy)LfL), (A.ll) 

where we combined the terms such that the order y is manifest. By the spectral theorem ( |A.11| ) 
is differentiable to any order in y, y > 0. Thus to complete the proof it suffices to show that 
the linear order of (|A.11|) vanishes. Since / — ~ ~y[H, /], we obtain 



(v) = ytT{ [H, f] {LfKfKfK - SKfLfKfK + 2LfLfKfK - SLfKfKfL 

-SKfLfLfK - 4LfKfLfK + LfLfLfK - SLfLfKfL)) + 0{y^). (A.12) 

We substitute L = K — 1 and use that [H, /]* = —[H, /]. The term with four K^s reads 

-mri[H,f]KfKfKfK) = 0. (A.13) 

The term with three K^s reads 

-6tT{[H, f]fKfKfK + KfKfK + KfKfK + KfKfKf)) = 0, (A.14) 

since first and fourth summand and second and third summand cancel. The term with two K^s 
reads 

-tr{[H, mfKfKf + 2KfK + 3fKfK + 3KfKf)) = 0, (A.15) 

since the first and second summand vanish and since the third summand cancels against the 
fourth one. The term with one K reads, upon adding the adjoint, 

tT{[H, f]f - f[H, f] - 3f[H, f]f + 3f[H, f]f). (A.16) 

At this point the specific structure of H enters. We have [H, f] = —f" — 2f\d/du). Working 
out the commutators ensures also that the last term vanishes. Thus the linear order from the 
4-cycles vanishes. The quadratic order does not vanish, however, and refiects the deviations 
from the local Brownian motion statistics. □ 

Appendix B: Trace class properties 

As shown in PaK is of trace class. If one establishes that R = e~y^ PbC^^ K is of trace class, 
then each summand in ( ^.9| ) is separately of trace class. In the energy representation R has the 
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kernel 



/oo 
dxki{x - A')Ai(x - A)e^^x(A < 0). (B.l) 

If A < 0, A' > 0, the integral is bounded uniformly in A, A'. If A < 0, A' < 0, one uses that 
Ai(x — A') < ce~i'^'~'^''^''^ for large negative A'. Thus the integral dominates the factor e~^'^ . 
This implies the bound 

|i?(A',A)| < ce-^(l^'l+l^l) (B.2) 
with suitable constants c, 7 > 0. Hence R is of trace class. 



Note added in proof: In their recent preprint p5[ Okounkov and Reshetikhin consider the 
(1, 1, 1) interface of the three-dimensional Ising model at zero temperature, which maps onto 
a domino tiling of a form similar to the Aztec diamond. They prove that correlations are of 
determinantal form and compute the limit shape. Their Theorem 1 is the analogue of our 
Eq. (|3§. 
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